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Abstract 

The study of non-stationary processes whose local form has controlled properties is 
a fruitful and important area of research, both in theory and applications. In [9], 
a particular way of constructing such processes was investigated, leading in particu- 
lar to multifractional multistable processes, which were built using sums over Poisson 
processes. We present here a different construction of these processes, based on the 
Ferguson - Klass - LePage series representation of stable processes. We consider vari- 
ous particular cases of interest, including multistable Levy motion, multistable reverse 
Ornstein-Uhlenbeck process, log-fractional multistable motion and linear multistable 
multifractional motion. We also show that the processes defined here have the same 
finite dimensional distributions as the corresponding processes constructed in [9]. Fi- 
nally, we display numerical experiments showing graphs of synthesized paths of such 
processes. 



Keywords: localisable processes, stable processes, Ferguson - Klass - LePage series rep- 
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1 Introduction 

This work deals with a general method for building stochastic processes for which certain 
aspects of the local form are prescribed. We will mainly be interested here in local Holder 
regularity and local intensity of jumps, but our construction allows in principle to control 
other properties that could be of interest. Our approach is in the same spirit as the one 
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proposed in |9], but it uses different metliods. In particular, in [9], multistable processes, 
tliat is localisable processes wliicli are locally a-stable, but where the index of stability a 
varies with time, were constructed using sums over Poisson processes. We present here 
an alternative construction of such processes, based on the Ferguson - Klass - LePage series 
representation of stable stochastic processes [TOl [131 [H] . This representation is a powerful tool 
for the study of various aspects of stable processes, see for instance [3l[T8]. A comprehensive 
reference for the properties of this representation that will be needed here is |19j . 

Stochastic processes where the local Holder regularity varies with a parameter t are in- 
teresting both from a theoretical and practical point of view. A well-known example is 
mult if r actional Brownian motion (mBm), where the Hurst index h of fractional Brownian 
motion [T2l IT5] is replaced by a functional parameter permitting the Holder exponent 
to vary in a prescribed manner [H El HH [16] . This allows in addition local regularity and 
long range dependence to be decoupled to give sample paths that are both highly irregular 
and highly correlated, a useful feature for instance in terrain or TCP traffic modeling. 

However, local regularity, as measured by the Holder exponent, is not the only local feature 
of a process that is useful in theory and applications. Jump characteristics also need to be 
accounted for, e.g. for studying processes with paths in -D(R) (the space of cadlag functions, 
i.e. functions which are continuous on the right and have left limits at all t ^T). This has 
applications for instance in the modeling of financial or medical data. Stable non-Gaussian 
processes yield relevant models for data containing discontinuities, with the stability index 
a controlling the distribution of jumps. 

Just for the same reason why it is interesting to consider stochastic processes whose local 
Holder exponent changes in a controlled manner, tractable models where the "jump intensity" 
a is allowed to vary in time are needed, for instance to obtain a more accurate description 
of some aspects of the local structure of functions in Z}(R). 

The approach described in this work allows in particular to construct processes where 
both h and a evolve in time in a prescribed way. Having two functional parameters allows 
to finely tune the local properties of these processes. This may prove useful to model two 
distinct aspects of financial risk, to describe epileptic episodes in EEG where for some periods 
there may be only small jumps and at other instants very large ones, or to study textured 
images where both Holder regularity and the distribution of discontinuities vary. 

Let us now recall the definition of a localisable process P, [Z]: ^ = {^(i) : t G R} is said 
to be /i— localisable at u if there exists an /i G R and a non-trivial limiting process such 
that 



(Note may and in general will vary with u.) When the limit exits, Y^ = {Y^{t) : t G R} 
is termed the local form or tangent process of F at u (see [21 [16] for similar notions). The 
limit (11. II) may be taken in several ways. In this work, we will only deal with the case 
where convergence occurs in finite dimensional distributions (equality in finite dimensional 

distributions is denoted ^=). When convergence takes place in distribution, the process is 

called strongly h-localisable (equality in distributions is denoted =). 

As mentioned above, a now classical example is mult if r actional Brownian motion Y which 
"looks like" mdex-h{u) fractional Brownian motion close to time u but where h{u) varies, 
that is 
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where Bh is index-h fractional Brownian motion. A generalization of mBm, where the Gaus- 
sian measure is replaced by an a-stable one, leads to multifractional stable processes, where 
the local form is an /i(M)-self-similar linear a-stable motion [201 l2T] . 

The h-local form at u, if it exists, must be /i-self-similar, that is Y^{rt) = r^Y^{t) 
for r > 0. In addition, as shown in [6l[7], under quite general conditions, Y^ must also have 
stationary increments at almost all u at which there is strong localisability. Thus, typical local 
forms are self-similar with stationary increments (sssi), that is r~^{Y^{u + rt) — Y^{u)) = Y^(t) 
for all u and r > 0. Conversely, all sssi processes are localisable. Classes of known sssi 
processes include fractional Brownian motion, linear fractional stable motion and a-stable 
Levy motion, see O [19]. 

Similarly to [9] , our method for constructing localisable processes is to make use of stochas- 
tic fields {X(t,v), (t,v) G R^}, where t is time, and where the process t X{t,v) is 
localisable for all v. This field will allow to control the local form of a 'diagonal' process 
Y = {X{t,t) : t G R}. For instance, in the case of mBm, X will be a field of fractional 
Brownian motions, i.e. X{t,v) = Bh{v){t), where h is a smooth function of v ranging in 
[a,b] C (0, 1). This is the approach that was used originally in [Ij for studying mBm. From 
a heuristic point of view, taking the diagonal of such a stochastic field constructs a new 
process with local form depending on t by piecing together known localisable processes. In 
other words, we shall use random fields {X{t, v) : [t, v) G R^} such that for each v the local 
form X'^{-, v) oi X{-,v) at V is the desired local form Y^ of Y at v. An easy situation is when, 
for each v, the process {X{t,v) : t G R} is sssi, since this automatically entails localisability. 

It is clear that, in this approach, the structure of X{-,v) for f in a neighbourhood of u 
will be crucial to determine the local behaviour of Y near u. A simple way to control this 
structure is to define the random field as an integral or sum of functions that depend on t 
and V with respect to a single underlying random measure so as to provide the necessary 
correlations. 

General criteria that guarantee the transference of localisability from the X{-,v) to Y = 
{X{t,t) : t G R} were obtained in [H]. We will make use of the following one: 

Theorem 1.1 Let U he an interval withu an interior point. Suppose that for some < h < t] 
the process {X{t,u),t G U} is h-localisahle at u E U with local form X'^{-,u) and 

P{\X{v,v) - X{v,u)\ > Iw-ul") -^0 (1.3) 

as V ^ u. Then Y = {X{t,t) : t G U} is h-localisahle at u with yj^(-) = 

In the sequel, we shall consider specific classes of random fields and use Theorem 11.11 to 
build localisable processes with interesting local properties. As a particular case, we will study 
multifractional multistable processes, where both the local Holder regularity and intensity of 
jumps will evolve in a controlled manner. 

The remaining of this article is organized as follows: we first collect some notations in 
section [2l We then build localisable processes using a series representation that yields the 
necessary flexibility required for our purpose. We need to distinguish between the situations 
where the underlying space is finite (section [3]), or merely a— finite (section Hj). In each case, 
we define a random field depending on a "kernel" /, and give conditions on / ensuring local- 
isability of the diagonal process. We then consider in section O some examples: multistable 
Levy motion, multistable reverse Ornstein-Uhlenbeck process, log-fractional multistable mo- 
tion and linear multistable multifractional motion. Section [6] is devoted to computing the 
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finite dimensional distributions of our processes, and proving tliat tliey are tlie same as tlie 
ones of tlie corresponding processes constructed in [9j . Finally, section [7] displays graphs of 
certain localisable processes of interest, in particular mult if r actional multistable ones. 

Before we proceed, we note that constructing localisable processes using a stochastic field 
composed of sssi processes is obviously not the only approach that one can think of. It is for 
instance possible to follow a rather different path and construct localisable processes from 
moving average ones by imposing conditions on the kernel defining the moving average. See 
[H] for details. 



2 Notations 

We refer the reader to the first chapters of |l9j for basic notions on stable random variables 
and stable processes. In particular, recall that a process {X{t) : t G T}, where T is generally 
a subinterval of R, is called a-stable (0 < a < 2) if all its finite-dimensional distributions 
are a-stable. Many stable processes admit a stochastic integral representation as follows. 
Write S'a((T, /?, /i) for the a-stable distribution with scale parameter a, skewness (3 and shift- 
parameter fi] we will assume throughout that fi = 0. Let [E, S, m) be a sigma-finite measure 
space. Taking m as the control measure and /? : — > [— 1, 1] a measurable function, this 
defines an a-stable random measure M on £^ such that for A E S we have that M{A) ~ 

(m(A)i/",/^/3(a;)m(c/x)/m(A),0). If /3 = 0, the process is termed symmetric a-stable, or 
SaS. 

Let 

J^a = ^a{E,£,m) = {/ : / is measurable and ||/||a < oc}, 
where || ||a is the quasinorm (or norm if 1 < a < 2) given by 



(/^|/(x)rm(rf^))^/" (a^l) 
\f{x)\m{dx) + \f{x)(3{x) In \f{x)\\m{dx) (a = 1) 



(2.4) 



The stochastic integral of / G J-'a{E,£,m) with respect to M then exists [19, Chapter 3] 
with 



/(/)= / fix)Midx)r^S^iaf,l3f,0), (2.5) 
Je 

where 

llfll ^ //(x)< ">/?(x)m(rfx) 



and a<^> = sign(a)|a|^, see [19], Section 3.4]. In particular, 

^\j(f)w^l4oi,P,p)\\fK {0<p<a) , 
' ^•^ \ oo (p > a) 

where c{a,P,p) < oo, see [IHl Property L2.17]. 

In this work, we will only consider the symmetric case, i. e. we take (3 = for the remaining 
of the article. We believe most results should have a counterpart in the non-symmetric case, 
although the proofs would probably have to be much more involved. 
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3 A Ferguson - Klass - LePage series representation of 
localisable processes in the finite measure space case 

A well-known representation of stable random variables is the Ferguson - Klass - LePage 
series one [H [TOl [131 El [IB]- This representation is particularly adapted for our purpose 
since, as we shall see, it allows for easy generalization to the case of varying a. 
In this work, we will use the following version: 

Theorem 3.2 ('^ Theorem 3.10.1]) 

Let {E,S,m) be a finite measure space, and M be a symmetric a-stable random measure 
with a e (0, 2) and finite control measure m. Let (rj)j>i be a sequence of arrival times of a 
Poisson process with unit arrival time, (l^i)j>i be a sequence of i.i.d. random variables with 
distribution m = m/m{E) on E, and (7i)i>i be a sequence of i.i.d. random variables with 
distribution Pi^ji = 1) = P{'ji = —1) = 1/2. Assume finally that the three sequences (rj)j>i, 
{Vi)i>i, o^nd (7j)i>i are independent. Then, for any f G J-'a{E , S , m) , 

„ oo 

/ f{x)M{dx) ^ {C^m{E)f"Y.^,T;'/y{V,), (3.7) 
Je 

where Ca '■= r(2-a)^cos(7ra/2) a 7^ 1, Ci = 2/7r (Theorem 3.10.1 in [19j is more general, 
as it extends to non-symmetric stable processes, that are not considered here.) As mentioned 
above, a relevant feature of this representation for us is that the distributions of all random 
variables appearing in the sum are independent of a. We will use (13. 7p to contruct processes 
with varying a as described in the following Theorem. 

Theorem 3.3 Let {E,S,m) be a finite measure space. Let a be a function defined on 
R and ranging in [c,d\ C (0,2). Let h be a function defined on R. Let f(t,u,.) be a 
family of functions such that, for all (t,u) G R^, f(t,u,.) G J-'a(u){E,S,m) . Let (rj)j>i be a 
sequence of arrival times of a Poisson process with unit arrival time, {Vi)i>i be a sequence 
of i.i.d. random variables with distribution m = m/m{E) on E, and (7j)i>i be a sequence 
of i.i.d. random variables with distribution P{'~fi = 1) = P(7i = —1) = 1/2. Assume finally 
that the three sequences {Ti)i>i, (Vi)j>i, and (7j)j>i are independent. Consider the following 
random field: 

oo 

X{t, u) = 6(n)(m(i?)) ^ 7.r-^/"^"V(t, n, V^, (3.8) 

i=l 

where Ca = [J^ x'"" sm{x)dx) \ Assume that X{t,u) (as a process int) is localisable at u 
with exponent h G (0, 1) and local form X'^{t,u). Assume in addition that: 

• (CI) The family of functions v — > f{t, v, x) is differentiable for all {v, t) in a neighbour- 
hood of u and almost all x in E. The derivatives of f with respect to v are denoted by 
f ■ 

• (C2) There exists e > such that: 

sup / sup (|/(t,«;,x)|"(^)) m{dx) < oo. (3.9) 

t£B{u,e) J E wGB{u,£) 
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• (C3) There exists £ > such that: 

sup I sup (|/^(t,w,x)|"("'))m(dx) < oo. (3.10) 

t&B(u,£) J E w&B{u,e) 



(C4) There exists £ > such that: 



sup / sup 

teB{u,e) J E weB{u,£) 



\f{t,w,x) log \f{t,w,x)\\'''^'"A m{dx) < oo. (3.11) 



Then Y{t) = X{t,t) is localisable at u with exponent h and local form Y^{t) = X'^{t,u). 
Proof 

The function u ^ C^a[u)^^ since a{u) ranges in [c,d] C (0,2). We shall denote 

a{u) = 6('u)(m(i?))-'^/"*^")C^^"|"''. The function a is thus also C^. We want to apply Theorem 
1.1. With that in view, we estimate, for v e B{u,s) (the ball centered at u with radius e), 

oo oo 

X{v, v) - X{v, u) =: ^iiUv) - M^)) + E ^^(^^(^) - ^^(^))' 

i=l i=l 

where 

$,(^) = a(ii)rV"H/(^,x,,y,) 

and 

The reason for introducing the $i and the is that the random variables Fj are not in- 
dependent, which complicates their study. We shall decompose the sum involving the 
into series of independent random variables which will be dealt with using the three series 
theorem. The sum involving the ^'j will be studied by taking advantage of the fact that, for 
large enough i, each Fj is "close" to i in some sense. 

In the sequel, since only the values of a inside B{u, e) matter, we shall agree by conven- 
tion that c denotes in fact inf^,gB(„ a(v) and likewise d — sup^g5(„^£) q;(v). Note that, by 
decreasing e, d — c may be made arbitrarily small. 

Thanks to the assumptions on a and /, $i and ^'^ are different iable and one computes: 

^l{u) = a'(«)i-^/"(")/(^, u, Vi) + aiu)i-'/''^^^f:,{v, u, V^} + a(«)^ log(i)i-V"(«)/(^, u, V^, 
and 

= a\u) (FrV«H _ rV"W) /(^, + a{u) (f-^/'^^") - r^-H) /^(^, 

(log(F,)F-^/"(") -log(i)rV«H) f{v,u,V^. 

The mean value theorem yields that there exists a sequence of independent random num- 
bers Wi G [u,v] (or [v,u]) and a sequence of random numbers Xi G [u,v] (or [v,u]) such 
that: 

oo oo 

X{v, u) - X{v, v)^{u- v) + Z] + Zl) + {u-v) Y,{Yl + + F/), 

i=l 2=1 



where 



a{Wiy 



Y- =-iia{x,, 

a[XiY 

Note that each Wi depends on a, f,a,u,v,Vi, but not on 7j. This remark will be useful in 
the sequel. 

The remainder of the proof is divided into four steps. The first step will apply the three- 

oo 

series theorem to show that each series Yl^i^j ~ 1)2,3, converges almost surely. In the 

i=l 

oo 

second step, we will prove that Yl ^/ ^-Iso converges almost surely for j = 1, 2, 3. In the third 

i=l 

oo 

step we will prove that condition (11.31) is verified by ^ Z^,j = 1, 2, 3. Finally, step four will 



i=l 



prove the same thing for X] — 1; 2, 3. 



1=1 



First step: almost sure convergence of '^Z!!,j = 1,2,3. 

i=l 

oo 

Consider Z^ = ^ Z}. Fix A > 0. We shall deal successively with the three series involved 

i=l 

the three-series theorem. 

oo 

First series : Si = Y. ^{\Z}\ > A). 

i=l 



< P| |/(i., tt..,V-.)|»>""' > ! inf 

w£B(u,£) 



A 



\a'[w) 



a{w) 



Note that, since a' is bounded on the compact interval [u,v], K := inf^ugB(„_£) 
is strictly positive. 



W{w)\ 



P(|Z^|>A) < P sup \f{v,w,V,)\''^'"^ > Ki 

\w£B(u,e) 

= P( sup \f{v,w,Vi)\''^'"^>Ki 
\weB{u,£) 



Thus 



+00 



+00 



^P(|Z/|>A) < 5^P( sup \f{v,w,V,r'-^^ > Kt 

1 



1=1 



sup i/(t.,tt.,Vi)r'"'' 

«)S-B(w,e) 



< 



sup / sup 

teB{u,£) J E w&B{u,e) 



{\f{t,w,x)r^'"^)m{dx) 



< +00 



by (C2). 

n 

Second series : = E E(^/l{|^/| < A}). 

E{Zll{\Zl\ < X}) = Ei^,a\w.,)t-'/''^'"^^fiv,w.,,V,)l{\a\w,)t''^''^^^^ 

= E(7i)E(a'K)^"'/"^"'V(^,^i, V^i)l{|a'K)ri/"('"')/(t;,«;„\/,)| < A}) 
= 0, 

where we have used the facts that 7^ is independent of Wi, Vi and E(7j) = 0. As a 
consequence, hm„^+oo S2 = 0. 

00 

Third series : The final series we need to consider is 5*3 = ^ E[(Z/1{|Z/| < 1})^]. Take 
A =10. 

Let f] be such that d < r] < 2. 



izli{\zl\<i}r < \zl\n{\zl\<i} 

E[{Zll{\Zl\<l}f] < E[|^/ri{|^^|<l}] 



+00 



p{\zl\n{\zl\<i}>x)dx 



= [\{\zl\n{\zl\<i}>x)dx 

Jo 

< [ P{\Zl\'' > x)dx. 
Jo 



^Recall that, in the three series theorera, for the series X^i^i -^i to converge almost surely, it is necessary 
that, for all A > 0, the three series J^Zi > >')^T,Zi E(X,l(|^i| < A)), and E^i Var(Xa(|^i| < A)) 

converge, and it is sufficient that they converge for one A > 0, see, e.g. [17], Theorem 6.1. 
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Now, for all x in (0, 1), 

P{\Zl\^>x) = P{\Zl\>x'/^) 



X "I 



< P( sup \f{v,w,Vi)\'^^'"^ > K'ix'^ 
P ( sup > K'ixn ] , 

\weB{u,s) 



where K' := inf^gB(«,£) 



|a'(«;)| 



is strictly positive. Thus, 



Jo y«)es(w,£) y 

< r^E( sup 

M sup / sup (|/(t,^,a;)r(-))m(cia;)Ur^ 

^ \teB{u,e) J E weB{u,e) J \J Xv 



< 

- K 
< +00 



The case of the Z'^ — ^ Z? is treated similarly, since the conditions required on (a', /) in 



i=l 



the proof above are also satisfied by (a, /^). 
Consider finally = ^ Zf. Fix A > 0. 



i=l 



First series : = £ Pd^f] > A). 



i=l 



p(|zf|>A) = p |/K^„y,)|> 



Aa(wi)2ii/°("'') 
|a(wj)Q!'(wi)| logi 



< p(|/K«;,,y,)r("'^)>x" 



(log i)oc{wi) J ' 



where := inf^eB(„,e) 

a'. In the sequel, K wil 
from line to line. 



a{w)a'{w)\ I strictly positive by the assumptions on a, a and 

always denote a finite positive constant, that may however change 
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Let gi{x) = (-iQg^-)a(z^i) for X > 1 and i eW. For x large enough and for all i, Qi is strictly 
increasing and limx-^+oo Qiix) — +00. For z large enough (independently of i), 



^^^^ ^ ^ ' {\ogz + a{wi) loglog;2)«K) 

Kz 



> 
- 2 



Let A > e be such that: ^z > A^^i e N\g-^{z) < f (log 
LetUi^\f{v,Wi,Vi)\-(^^\ 



P(|Zf I > A) < P{Ui > K- — ^—-^) 
^' " ' ^ (logz)"("'')^ 



< P( [{u, <A}n {u, > K^^^^}^ u ({[/, > .1} n {g-\u;) > i})) 

< P{A >Ui> K-^^) + P{{Ui > A} n {Ui\ logC/,r("'^) > Ki}) 

< P{A>Ui>K—^) + P{U^\\ogUi\''^^^^>Ki) 

[logiY 



< P{A>Ui>K——^) + P{ sup [\f{v,w,V,)\og\f{v,w,V^)r^'"^]>-i). 
On the one hand, P{A > Ui > K^^^) vanishes for i large. On the other hand, by (C3), 

^P( sup [\f{v,w,V,)log\f{v,w,V,W^'"^]> ^t] < 



^E( sup [\f{v,w,V,)log\f{v,w,V^W^^^]\ <+oo 



K 



and thus < +00. 

n 

Second series : 5^ = E E(^fl{|^f | < A}). 

i=l 

For the same reason as in the case of Z] [i.e. ji is independent of Wi, V^), 5*2 = 0, Vn. 
00 

Third series : S3 = E{{Zfl{\Zf \ < l})^). Take again A = L 

i=l 

Let ?7 > be such that d < i] < 2, and > be such that < //. < 1 — ^. Using the same 
line of reasoning as in the case of Z^, 



E(Zf I < l}f < f P{\Zf\^ > x)dx. 

Jo 
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We have, for all x in (0, 1), 

P{\Zf\'^>x) = P{\Zf\>x'^'^) 

< p(\f{v,Wi,Vi}\''^^'^>Kxv 



where K := mi^^B{u,e) 
a'. 



\a(w)a' {w)\ 



is strictly positive by the assumptions on a, a and 



P{\Zf\^>x) < P{A>Ui>Kxl-^) + P{ sup [\f{v,w,Vi)\og\f{v,w,Vi)\r^^^]>^'' 



(logi) 



Let us deal with the first term of the sum in the right hand side of the above inequality. 



P{A > Ui > Kxl —^—;) < P{A> Ui> Kxvi^-") 



{logiY 



for i large enough, i.e. i > i*, where i* depends only on d and /i. As a consequence : 



+00 



P{A >Ui> Kx^ 



d I 



(logi) 



;\d' 



< 



< 



But 



< 1 and thus 



+00 



A \i 



d 

Kx^ 



K' 



fi +00 

is 



P(A >Ui> Kxv- —)dx < +00. 



Now for the second term in the sum: 

+00 



d 

Kx-^ 



Y,P{ sup [\j^v,w,V^)\og\f{v,w,V^W^'"\>—-i) 



< J!_E( sup [\f{v,w,V,)\og\J{v,w,VM^^'"'^] 
Kx"^ \weB{u,£) 



< d 



K' 

d 

X'n 
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and as a consequence 

/ ^P( sup [\f{v,w,V^)\og\f{v,w,V,W^'"^]>-—^)dx 

Jo weB(u,£) " 



K' , 

< / —dx 

Jo 

< +00. 



We have thus verified all the conditions in the three series theorem, and shown that the 
series Z^, and are almost surely convergent. 

00 

Second step: almost sure convergence of Yl^i^j ~ 1)2,3. 

i=l 



To prove that the series Yl^i^J — 1)2,3 converge almost surely, we will first show that 

i=l 

it is enough to prove that X] ^/l{i<£i<2}n{|y.^|<i} converges almost surely for j = 1,2,3. 

1=1 2— I — I — 

Indeed, we prove now that Ylili ^ {{j < x — 2} U {|^/| > l}j < 00, where T denotes the 



complementary set of the set T. 



P({^<y <2}U{|r/|>l}j = p({i<^<2}U 



{|>^^l>i}n{^<^<2} 



< P({r,: < ^}) + P({r. > 2^}) + p > 1} n < ^ 



Fj, as a sum of independent and identically distributed exponential random variables with 
mean 1, satisfy a Large Deviation Principle with rate function A* (a;) = x — 1 — log(x) for 
X > and infinity for x < (see for instance |1] p. 35), thus Yl ^({^i < f}) < +^ ^^"^ 

i>l 

E P({r, > 2^}) < +00. 



i>l 



Consider now ^ P > 1} n {| < ^f < 2}) , for j = 1, 2, 3. 



i>l 



Case j = 1 or j = 2 : 



p(^{|F/|>i}n{^<^<2}) = p({|^fll(y)-^/"^^'^)-i|>i}n{i<^<2}) 

< P |^{(2V"(-0 _ i)\zi\ > 1} n < ^ < 2}^ 



< P \z'\> 



21/d - 1 



thus E.>i P > 1} n {i < < 2}) < +00. 
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Case j — 3: 
For i > 1 : 

P({|^/I>i}n{^<^<2}) = p({|^f||i^(^)-^/"(-)-i|>i}n{^<^<2} 

~ l)(i>i) bounded on {| < ^ < 2}, thus there exists K >0 such that 

P({|>^/I>i}n{^<^<2}) <p(^|zf|>-i 

which entails that X) P ({l^i^l > 1} n {| < < 2}) < +oo. 

i>l * 



We are thus left with proving that 5^ 5^''l^i<ri^2}n{|y.^|<i} converges almost surely for 
J = 1,2,3. 

In that view, we shall apply the following well-known lemma: 
Lemma 3.4 Let {Xk,k > 1} be a sequence of random variables such that E|X„| < +00, 

n=l 

+00 

then ^ Xn converges almost surely. 

n=l 

00 r . -1 

Let us show that E l^''|l{i<£i<2}n{|Y.^|<i}J < +°°- 



< jy(^{\Yf\>x}n{\<^<2}yx. 



Let fi. = {i < ^ < 2}. 
Case J = 1: 



_ 1 -. 

Using the finite- increments formula applied to the function x ^ x on [3,2], one 
easily shows that 

P{{\Y^\>x}nB,) < p(^{\zl\\^-i\>x^JnB^ 

= P (^\a'{x,)i-'/-^^^^f{v,Xi,Vi)\\^ -1\> x-^JnB, 
< P (^{\f{v,Xi,Vi)\'^^''^^\j -1]''^''^^ > K^ix"^''^^} n B,^ 
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where := ini^^B{u,e) 
a. Thus 



2l+i/=|a'(«;)| 



a{w) 



is strictly positive by the assumptions on a' and 



P({|r/| >a;}n5,) < P( sup - 1|" > iTcix" 

Case d>l : 

Fix 7/ G ((i, 1 + |) (since a is continuous and c? < 2, by decreasing if necessary £, one may 
ensure that o? < 1 + c/2). By Markov and Holder inequalities, 



P({|F/| >a;}nS,) < P 



sup 



< 



< 



1 



K 1 







) 




r. 


c/2r, / 


E|--l|' 




i 





1/r; 



1/7? 



sup E( sup |/(^;,'i/;,T/i) 

v(iB{u,£) w€B(u,s) 



where we have used that the variance of Fj is equal to i, and K does not depend on v thanks 



to assumption (C2). Thus E 
Case d < 1: 



i^{^<ii^<2}n{iy/|<i} 



< -UT^ where ^ + |^ > 1- 



P({|F/| >x}nA) < ^i^E( sup \f(v,w,V^T^'"^)E\^-l\ 

X IS-cl weB{u,e) * 



< i^J_l(E|I^-l|2)c/2 



thus E 



'■^^^'■'■{^<^i<2}n{|y.i|<i} 



< 



jl+c/2 ^d' 

^ with 1 + S > 1. 



jl+c/2 



The case of ^ E 



\^^\'^{B,n{\Y-'\<i} 



is treated similarly, since the conditions required on 



(a', /) in the proof above are also satisfied by (a, /^). 
Case j — 3: 

We now consider ^ E \Y^\1^b.^^\y^\^i} ■ 

j>l L 'J 

Again by the finite-increments formula, there exists a constant i^^c.d) which depends on c 
and d, such that, for i > 1, 



log j Z 



i 



lB,<K,,d\--l\lB,. 



14 



Then, 



p({|r,=^| >x}n5,) < p(^{\zf\K,^,\^-i\>x}nB,^ 



< P| sup \f{v,w,Vi)\'^^'"^\— -11" > K- —x" 



wGB{u,£) 



(logi) 



Case d> 1 : 
Fixr/e (t^,l + f). 



P(m>x}nB-) < ^(!^ 



Case d < 1 



P{{\Yf\>.}nB.) < ^ 



d ^l+c/2 



X' 



+00 

As a conclusion, for j = 1, 2, 3, Yl ^1 converges almost surely. 

We now move to the last two steps of the proof: to verify /i-localisability, we need to 
check that for some r] such that h < t] < 1, P{\Z^ > \v — u\^~^) and Pdl^-'l > |f — ul"^'^) 
tend to when v tends to u, for j = 1, 2, 3. 

Third step: verification of ( 11.31) for Z^,j = 1,2,3. 

We need to estimate Pd^-'l > \v — u\^'^^). 
LetaG(0,^). 



oo 



\V — U\ 



i=l / \ i=l 

oo 



\Zi\>\v-u\- 
i=l 



> 



Since 7j is independent from 7^ for i 7^ A; and is independent of 7^, Markov inequality 
yields 



1=1 



Let 7 G {d,2). For any M > 1, using the same computations as in the first step, third 
series (page [TOl) . we get: 
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\Zl\<M 



■J 12 



< MM P(|Z/| > Mx^/'<)dx. 
Jo 



M2 l^'l<*^ 

For j = 1, 

[ P{\Zl\> Mx^/'')dx < [ P( sup \f{v,w,Vi)\''^'"^ > MHKx'^/^)d. 

Jo Jo w&B(u,£) 

thus 

oo 

1=1 

The same conclusion holds for j = 2: 

oo 



|Zi|<A/ 



oo „i 



Jo weB{u,e) 



i=l 



< KM 



2-c 



For j = 3, choose fi E (0, 1 — -), and compute as on page fTTl 



P{\Zf\ > Mx^''<)dx < I P{A> \f{v,Wi,Vi)\''^'"^'^ > KM''x'^'^—^)dx 



+ fp{ sup [\f{v,w,V,)\ogf{v,w,V,)\r^'^'^>K 

Jo w£B{u,e) 



veB{u,£) 

For i large enough, i > i* where i* depends only on d and /x, 



d'i ' 



thus 



log(i)' 



i=l 



log(z)' 



< 



E^l^ s^P E( sup 



,i=l 



v£B(u,e) w^B(u,e) 



As a consequence, for M > 1: 



]^ oo 



i=l 



5^PU>|/(t;,w„\/,)r(-')>irM 



log(0^ 



< KM-" + KM- — 



< KM-". 
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Since 

J2 P( sup [\f{v,w,V^)log\f{v,w,V,)\\r^^^>K—^i)dx<KM-'^, 

Jo weB{u,£) " 



we get 



oo 

1=1 



\Zl\<M 



Let M —\v — u\ "■. Using previously obtained inequalities, we get, for j — 1, 2, 3: 

/ oo I _ |r;-l\ 

P I I V ■ I > 1^ "I < K\v - 7/|2(l-'?)-«(2-c) 



and 

lim P ( I Vz/l|7,,,, i-al 1^0. 



\V — u\'^ ^ 



We 



2 

i=l 



consider now the second term P ^| Yl^i ^i-^\z^\>\v-u\-<^\ 2 
Let i* = mi{n > 1 : i > n, < \v - Since ^ P(|Z/| > - < +00, the 



Borel-Cantelli lemma yields P{i* = +00) =0. As a consequence. 



Iri— 1 \ 00 / 00 I In — 1 

\V — U\ ' \ ST^ ^ I r, „i , — U\ 



= E P I {l E ^Aiz^l>|.-.|-^l > ^^^^} n {^* = n} 

n=2 \ 1=1 

00 



i=l / n=l \ i=l 

00/00 



n=2 



For n > 2, P(i* = n) < P(|^^_i| > It; - 

For j = 1, P(i* = n) < P(sup^gB(„^^) > |-t; - u\-^^K{n - 1)), and thus 



J2P{i* = n) < Klv-ul'^'Ei sup |/(^;, w, 

n=2 weB{u,e) 

< Klv-ul"" sup E( sup 
teB(w,£) iueB(u,£) 



For j = 2, 



^ P(i* = n) < - sup E( sup \fl{t,w,Vi)\''^'^^] 

„_2 teB{u,e) ■weB{u,e) 



17 



and for j = 3, 



Y^P{e = n) < K\v-ur sup E( sup \f{t,w,Vi)\og\f{t,w,V,W^^^] 

t&B{u,£) w€B{u,e) 

n>2 



\og{i 



We have shown previously that the second term in the sum on the right hand side of the 
above inequahty is bounded from above hy K\v — u\°'^ . Finally, 



W — u\ 



0. 



i=l 



Fourth step: verification of ( 11.31) for Y^,j = 1,2,3. 
We consider now PdV-'l > |f — u\''^^^). 

Let i* = inf{n > 1 : i > n, \Y/\ < 1 and ^ < ^ < 2}. Since ^ Pi{\Y/\ > 1} U {F^ < 

* j>i 

|} U {Fj > 2i}) < +00, the Borel-CanteUi lemma yields P{i* = +oo) =0. As a consequence, 

n>l \ i=l / 

Let bn{v) = P ^{| X] ^/l — 1"^ n {i* = n}^ . Our strategy is the following: we show 

that, for each fixed n, bn{v) tend to when v tends to u. Then we prove that there exists a 
summable sequence (c„)„ such that, for all n and all v, bn{v) < c„. We conclude using the 
dominated convergence theorem that J2n>i ^n{v) tends to when v tends to u. 
For all n > 1, 



i=l 



n—l 

For n>2, consider P({| E ^/l > f ))• 

1=1 



i=l i=l 



2(n- 1) 



Let p G (0, ^). With if a positive constant that may change from line to line and depend 
on n but not on v, we have, for j = 1: 
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pI\YI\> ^1 ^' < P sup |/K^,T/,)r(-)|(i^)-V'^(-')-l|«(-^)> ^'"^ 



(2(n - l)a'(a;i))°(^^) 



< P I ( sup V^i)|"('"))f|(— )-^/"("') - > K\v - u\P'^''-^'> 

weB(u,£) ^ 



Using the independence of Vi and Fj and Markov inequality, 

p (\Yl\ > < i^|^;-«|f'=a-'?)E ( ( sup \f{v,w,Vi)\''^''^y] E(|(— - ir(^^)) 

\ 2(n - 1) y y «,eB(w,£) J I 

r . 

~ i 

It remains to check that the expectation in the right hand side of the above inequahty is 
finite: 

F F F 

X % X 

< 1 + E(|(— )-^/"(^') - l|*'"(^')lr,<i) 

< l + E(|(^i^)-^/^-ir(-')lr,<0 

< 1 + E(i(— - in + E(i(— )-^/" - 1|^'^). 

Since p < % E(|(^)-i/^ - Ip^) < +00 and E(|(^i)-^/'= - l]?"^) < +00 (this is easily verified 
by computing these expectations using the density of Vi). Thus we have E(|(?^)~^/"(^') — 
< +00, and 

\v — u\'^~^ 



hmP |y/|> ' , HO- 

Since the conditions required on (a', /) are also satisfied by (a, /^), 



'7-1 ■ 

lim P ( IK^I > '\ . 1 = 0. 



.21- l-f-^l 



2(n- 1) 

We consider now the case j = 3. When i — 

Nl^.'l^^) - P (1 iog(r.)rr'^°'-'/(.. X., v-OI > - 1) - 



(iog(ri))"(^i) ^^ 
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where again K depends on n but not on v. Since p < 1 and a is bounded, E ^( 
+00, and 



limP(|F3|>tL_JitL! 
v-^u V ' ^ ' - 2(n - 1) 



0. 



For i > 2, 



P |F/| > 



\v — u 



r,-l 



2(n- 1) 



> 



\og{{)2\a{xi)a'{xi)\{n - 1) 



< ^ , I _ Apc{l-v)n ( | l0g(r») .rj^-i/a(a;,;) _ |a(a;i)p 



/ l0g(ri) Ti _ ^ 



iog(r,) T, 



,iog(r,) T, 



log(i) ^ i ^ ' ^ ' log(i) ^ i > ' 



^1 |iog(r\)^ri^_,/,_^l„(^^) \ ^ Ef|i^^^(^^)-^/'=-i|^f^ +Ef|i^^^(^^)-^/'=-i|* 



log(i) i 



log(i) i 



log(i) i 



and 



^£g(£0 /£i _ 1 \a{xi)p 

log(i) ^ i ^ ' 



log(i) i 



log(i) i 



Since p e (0> ^^i® four terms in the right hand sides of the two last inequahties are 
finite (use again the density of Vi) and thus E (\^-^^{^y^''^^'''^ - < +00. As a 

1,7-1 ■ 



consequence, 



hmP IK'' I > 



W — u\ 



2(n- 1) 



0. 



Finally, we have, for j e {1, 2, 3}, 



n-l 



limP({|^l^^|> })^0 



r,-l 



/ 00 

Let us now consider, for n > 1, P M | X) ^i'^{\Y^\<\}r^{\<U<2\\ - 



\v—u\ 



V-1 



\v — ul"^ ^ 



{|^'l<i}n{|<^f<2}l - 2 



} < 2|^;-m|^-''E 



{\Yi\<i}nO^<^<2}\ 



{|^^l<i}n{i<if<2} 



i=l 



< K\v-u\^-'' 
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(recall that the constants K used in bounding the series ^{\Yi\'^j^Y/^<^<i}n{^<^<2}^ '^'^^ 
depend on v). Thus hn{v) when f — > m for each n. 

In view of using the dominated convergence theorem, we compute (recall that -Bj = {| < 

1^<2}): 

hn{v) < P{{e = n}) 

< P{{\YU\>l}VJB:r,) 

< > 1} n + P(^ <l) + P(^ > 2). 



For j = 1 and d > 1, 



p({i>;M>i}ni?„-i)<^ — E( sup \fit,w,v,r^-^)y/^ 

and if d < 1, 

P({|>;-il>l}nS„-i)<7 ^^( sup E( sup \f{t,w,V,r^-^)). 

The same conclusion holds for j = 2, while, for j = 3, 

^'^ .(log(n - l))^/'?. 



This finishes the proof. ■ 



Ld>l- 



4 A Ferguson - Klass - LePage series representation 
of localisable processes in the cr-finite measure space 
case 

When the space E has infinite measure, one cannot use the representation above, since 
it is no longer possible to renormalize by m{E). This is a major drawback, since typical 
applications we have in mind deal with processes defined on the real line, i.e. E = and 
m is the Lebesgue measure. However, in the cx-finite case, one may always perform a change 
of measure that allows to reduce to the finite explained in |19j, proposition 3.11.3 

(for specific examples of changes of measure, see section [5]). In terms of localisability, this 
merely translates into adding a natural condition involving both the kernel and the change 
of measure: 



Theorem 4.5 Let {E,S,m) be a a-finite measure space. Let r : E ^ IR+ he such that 
rh{dx) = :^m{dx) is a probability measure. Let a be a function defined on R and 
ranging in [c,d] C (0,2). Let h be a function defined on R. Let f{t,u, .) be a family of 
functions such that, for all (t,u) G R^, f(t,u, .) G J-'a{u){E,S,m). Let (ri)j>i be a sequence of 
arrival times of a Poisson process with unit arrival time, {Vi)i>i be a sequence ofi.i.d. random 
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variables with distribution m on E, and (7i)i>i be a sequence of i.i.d. random variables with 
distribution P{ji = 1) = Piji = —1) = 1/2. Assume finally that the three sequences (Ti)i^i, 
{Vi)i>i, and {'~fi)i>i are independent. Consider the following random field: 

oo 

X(t,«) = 6HC^(3-)5^7.rr^/"("V(\/.)V"H/(t,«,y,), (4.12) 

1=1 

where Ca = (/q°° sin(x)(ix) ^ . Assume that X{t,u) (as a process in t) is localisable at u 
with exponent h e (0, 1) and local form X'^{t,u). Assume in addition that: 

• (Csl) The family of functions v — f{t,v,x) is differentiable for all {v,t) in a neigh- 
bourhood of u and almost all x in E. The derivatives of f with respect to v are denoted 
f. 

J V 

• (Cs2) There exists e > such that: 



sup 



sup (|/(t,ti;,x)r("'))m(rfx) < oo. (4.13) 



teB{u,e) J E w£B(u,e) 

(Cs3) There exists e > such that: 



t£B{u,£) J E w&B{u,£) 

( Cs4 ) There exists e > such that: 



sup [ sup {\f^{t,w,x)\''^'^^) m{dx) <oo. (4.14) 

\B(u,£) J E w&B(u,e) 



sup / sup \f {t,w,x) log \f{t,w,x)\\°'^^^ m{dx) < OO. (4-15) 

t£B{u,e) J E w£B{u,6) 



(Cs5) There exists e > such that : 



p / sup |/(t, w, x) log(r(x))|"*-"'^ m{dx) < oo. (4.16) 

u.e) J E w£B(u,e) '- ■' 



sup 

t&B{u,e) J E weB{u,£) 

Then Y{t) = X(t,t) is localisable at u with exponent h and local form Yl^{t) = X'^{t,u). 

Remark: from (I4.12p . it may seem as though the process Y depends on the particular 
change of measure used, i.e. the choice of a specific r. However, this is not case. More 
precisely, proposition 16.111 below shows that the finite dimensional distributions of Y only 
depend on m. 

Proof 

We shall apply Theorem 13.31 to the function g{t, w,x) = r(x)i/"("')/(^, 

w, x) on {E, S, rh). 

• By (Csl), the family of functions v f{t,v,x) is differentiable for all {v,t) in a 
neighbourhood of u and almost all x in E thus v — > g(t,v,x) is differentiable too and 
(CI) holds. 
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Choose e > such that (Cs2) holds. 



sup 



One has 



sup m(rfx) 

R weB(u,e) 



w£B{u,e) 



r{x) sup (|/(t,w,x)|"('"))m(rfx) 
sup {\f{t,w,x)\''^'"^)m{dx) 

R ■weB{u,e) 



thus 

sup / sup {\g{t,w,x)\"^'^'') rh{dx) 

t£B{u,e) Jr weB{u,e) 



sup / sup {\f{t,w,x)\"'-'^'')m{dx) 

t£B{u,e) Jk weB{u,e) 



and (C2) holds. 

Choose e > such that (Cs4) and (Cs5) hold. 



sup 

R w£B{u,e) 



\g{t,w,x) \og\g{t,w,x] 



< / r(x) sup 

'r w£B{u,e) L 



< / sup 

Jr weB{u,£) 



\f{t, w, x) log |r(x)^/"('^)/(t, w, x) I f ^'"^j m(dx) 
/(t, tt;, x) log |r(x)^/"("'V(^, ^i^, f ^"'^1 m(da;). 



Expanding the logarithm above and using the inequality |a + < max(l, 2^ ^)(|a|'' + 
valid for all real numbers a, b and all positive 6, one sees that (C4) holds. 

Choose e > such that (Cs3) and (Cs5) hold. 

giit, w, x) = r(x)i/"(-) w, x) - ^ log(r(x))/(t, x] 



and 



< / sup 

'r «>e-B(u,e) 



sup {\g'^{t,w,x)r'^'^y)m{dx) 

R w£B(u,e) 



f^{t, w, x) ^— log(r(x))/(t, w, x] 



m{dx). 



The inequality \a + < max(l, 2^-^){\a\^ + shows that (C3) holds. 
Theorem 13.31 allows to conclude. 
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5 Examples of localisable processes 



In this section, we apply the results above and obtain some localisable processes of inter- 
est. In particular, we consider "multistable versions" of several classical processes. Similar 
multistable extensions were considered in [9], to which the interested reader might refer for 
comparison. 

We first recall some definitions. In the sequel, M will denote a symmetric a-stable 
(0 < a < 2) random measure on R with control measure Lebesgue measure C We will write 

L„(t) := f M{dz) 
Jo 

for a-stable Levy motion. 

The log-fractional stable motion is defined as 

/oo 
(log(|t-a;|)-log(|a;|))M(da;) (t G R). 
-oo 

This process is well-defined only for a G (1,2] (the integrand does not belong to J-'a for a < 
1). Both Levy motion and log-fractional stable motion are 1/a-self-similar with stationary 
increments. 

The following process is called linear fractional a-stable motion: 

/oo 
fa,H{b^,b~',t,x)M{dx) 
-oo 

where t G M, G (0, 1), 6+, 6" G M, and 

f^Mb^, b-,t, x) = 6+ ((t - x) J-^/" - (-x)?-^/") 

+ r ((t - x)^"^/" - (-x)^-^/") . 

La,H,b+,b- is again an sssi process. When 6+ = 6^ = 1, this process is called well-balanced 
linear fractional a-stable motion and denoted La,H- 
Finally, for A > 0, the stationary process 

/oo 
exp(-A(x - t))M{dx) [t G R) 

is called reverse Ornstein-Uhlenbeck process. 

The localisability of Levy motion, log-fractional stable motion and linear fractional a- 
stable motion simply stems from the fact that they are sssi. The localisability of the reverse 
Ornstein-Uhlenbeck process is proved in |8j. 

We will now define multistable versions of these processes. 

For the multistable Levy motion, we give two versions: one is fitted to the case where the 
time parameter varies in a compact interval [0,T], and one where it spans R. 

Theorem 5.6 (symmetric multistable Levy motion, compact case). Let a : [0,T] — [c, C 
(1,2) and b : [0,T] R'*" be continuously differentiable. Let (rj)j>i be a sequence of arrival 
times of a Poisson process with unit arrival time, (Vi)i>i be a sequence of i.i.d. random 
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variables with distribution rh{dx) , the uniform distribution on [0,T], and (7i)i>i be a sequence 
of i.i.d. random variables with distribution P{'~fi = 1) = P{'~fi = —1) = 1/2. Assume finally 
that the three sequences (ri)j>i, (Vi)j>i, and (7i)j>i are independent and define 

+ 00 

Y{t) = 6(t)Ci(;;«ri/"W 5^7.rr'/"^*^l[o,*](V^.) it e [0,T]). (5.17) 

1=1 

then Y is 1/ a{u)-localisable at any u G (0,T), with local form = h{u)La[u)- 

The proof is a simple application of Theorem 13.31 and is omitted. 

Theorem 5.7 (symmetric multistable Levy motion, non-compact case). Let a : R ^ [c, rf] C 
(1,2) and 6 : R — ^ R^ be continuously differentiable. Let (rj)j>i be a sequence of arrival 
times of a Poisson process with unit arrival time, (Vi)i>i be a sequence of i.i.d. random 
variables with distribution m{dx) = 2~-'ly_ij[(x)(ix on R, and {%)i>i be a sequence of 
i.i.d. random variables with distribution P{'~fi = 1) = -P(7j = —1) = 1/2. Assume finally that 
the three sequences (rj)j>i, (Vi)j>i, and (7j)i>i are independent and define 

yit) = EE7^rr^/"W2^-/aWij^_^j^^^._^^^.^(^^) e (5.18) 

i=l j=l 

then Y is 1/ a{u)-localisable at any u G R+, with local form Y^ = b{u)La{u)- 
Proof 

We apply Theorem 14.51 with m{dx) = dx, r{x) = Xljli 2''l[i-i,j[(a;), f(t,u,x) = l[o,t](x) 
and the random field 

oo 

Xit,u) = b{u)Cl{:^''^ Yl 7.r-^/-(")2^V.Hip^^j^j^_^_^.j(y^). 

X{.,u) is the symmetrical a('u)-Levy motion [19] and is thus ^^^-localisable with local 
foTmXi{.,u)=X{.,u). 

• (Csl) The family of functions v — *• f{t,v,x) is differentiable for all {v,t) in a neigh- 
bourhood of u and almost all x in E. The derivatives of / with respect to u vanish. 

• (Cs2) 

\f{t,w,x)r^^^ = i[o,t]{x) 

thus 

[ sup (|/(t,w;,x)|"("')) rfx = t 

and (Cs2) holds. 

• (Cs3) = so (Cs3) holds. 
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(Cs4) f{t,w,x)\og\f{t,w,x)\ =0 so (Cs4) holds. 
(Cs5) 



+00 

\f{t,w,x)log{r{x)r^''^ = 5^j"(-)log(2)"(-)l[o,]ny-i,[(x) 

i=i 

+00 

< log(2)'^^j'^l[o,i]nb-i,i[(x) 



thus 



sup 

R «)GB(m,£) 



\f{t, w, x) log(r(x)) r^"') rfx < log(2)'' f 



and (Cs5) holds 



Theorem 5.8 (Log-fractional multistahle motion). Let a : R — [c, rf] C (1,2) and b : 
R — s> R+ 6e continuously differentiable. Let (rj)j>i he a sequence of arrival times of a 
Poisson process with unit arrival time, (Vj)j>i he a sequence of i.i.d. random variables with 
distribution rh{dx) = ■^'^'^^j^'^'^[-j,-j+i[u[j~i,j[{x)dx on R, and (7i)i>i be a sequence of 
i.i.d. random variables with distribution P{'yi = 1) = P{'~fi = —1) = 1/2. Assume finally that 
the three sequences {Ti)i>i, (Vi)j>i, and (7i)i>i are independent and define 

+00 +00 2/a(t) 

yit) = ^ EE^''rri/"W(log|t-l^,|-log|\/,|)^j2^^^^^ {t e R) 

i=l j=l 

(5.19) 

then Y is l/a{u)-localisable at any m G R, with = b{u)Aa(u)- 
Proof 

We apply Theorem 14 . 5 1 wit h m{dx) = dx, r(x) = ^ X]jlii^l[-j,-j+i[u[j-i,i[(3^)) f(t^ x) = 
log{\t — x\) — log(|a;|) and the random field 

X{t, u) = b{u)Cl{:l-^ 7.rr^/"^"^(log It - V,\ - log |^d)TTT7^j'/"^"^lH-,-.+i[ub-i.[(V^.). 



X{.,u) is the symmetrical a(M)-Log-fractional motion. It is ^^^-localisable with local 
form X;(.,m) = b{u)Aaiu) [9J. 



(Csl) The family of functions v — *■ f{t,v,x) is differentiable for all {v,t) in a neigh- 
bourhood of u and almost all x in E. The derivatives of / with respect to u vanish. 
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• (Cs2) Va > 1, 3Ka > such that \f{t,w,x)\''dx < Ka\t\ so 

< 

and 

/ sup (|/(t,w,x)r(-Va^ < {K, + Kd)\t\ 

JR weB{u,e) 

thus (Cs2) holds. 

• (Cs3) = so (Cs3) holds. 

• (Cs4) 

|/(i,^/;,x)log(|/(i,^,x)|)r(-) = |/(t,w,x)log(|/(t,«;,a;)|)r(-)l||iog(|/(t,.,.)|)|<i} 

+ w, x) log(|/(t, a;)|)|"("')l{|iog(|/(t,^,x)|)i>i} 
< x) r^'") + x) log(|/(i, x) |)r(-)l{|^(t,^,,)|>e} 

^, x) log(|/(i, w, 

We shall bound each of the three terms that are added up in the right hand side of the 
above inequality. For the first term, 

|/(t,w,a;)|°('") < sup \f{t,w,x)\"'-'"\ 
For the second term, fix > 0, e > such that Vx > e, |a;log(|a;|)| < ii'|a;p"'"^ 

\f{t, W, X) l0g(|/(i, W, X) I) r(-)l{|^(,,^,,)|>e} < K\f(t, W, X) \'^'+^\ 

For the third term, fix i^i < < and > > such that 

^{\f{t,w,x)\<^} < l]-oo,Xi|t|[(a;) + 'i^]K2\t\,K3\t\[{x) + l]K4\t\,+oo[{x), 

then 

\f{t,w,x)log{\f{t,w,x)\)r^^\f^,^^^,)\^^ < \f{t,w,x)\'^\log{\f{t,w,x)\)\%_^,Km[i^) 

+ W, X) n log(|/(t, W, x)\)\%K2\t\,Kmi^) 

+ \f{t, w, a;)n log(|/(t, w, x)\)\%Kiit\,+ooi{x). 



log(|t-a;|) -log(|a;|)| 



iog|i--ir(^) 

X 

log|l--||'+| 



log|l 



rr 
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The function x i— > l^l^^l log(|x|)|'^ is bounded for x < K With M denoting an upper 
bound of this function, one has 



/ \f{t,w,x)\'\\og{\f{t,w,x)\)\^dx < M{K^-K2)\t\. 

JK2\t\ 



With M = 1 — -, we obtain 

"+~ /■+°° |log|M|niog|log|ii|||'^ 



i>+oo r+oo 

/ \fit,w,x)r\logi\f{t,w,x)\)fdx < K'\t\ + \t\ 



(1 - uY 
< o{\t\), 



du 



where K' and K" are numbers verifying K" > 1 and K' > 0. 

For the same reason, /_^'*' \ f{t,w,x)\'^\ log{\f{t,w,x)\)\'^dx < K\t\. We conclude that 



sup / sup 

teB{u,e) Jr weB{u,e) 



\f{t,w,x)log{\f{t,w,x)\)r^-^ 



dx < oo. 



(Cs5) 



\f{t,w,x)logirix))r^^^ < K,\fit,w,x)log{-)r^^^ 

+ir2j]i/(i,^,x)r(-)(iog(j))'^iH,-,+i[ub-i,,[(x). 

For j large enough {j > j*), _j+i[u[j-ij[(a;) < Xsj^li-^ -j+i[ub-ij[(a:;)). 

Thus 



|/(i,^,x)log(r(x))r(-) < i^6|/(^,^,x)r(-) + ir7j](log(jO)'^|;^l[-,W+i[ub-iJ[(^)- 



To conclude, note that 

dx 



[ (log(j))'^T^lH-,-,+i[ub-i,.[(^)^^ = 2(log(j))'^ r 

(log(j)) 



2(c-l)- 



c 
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Theorem 5.9 (Linear multistahle multifractional motion). Let 6 : R ^ , a : H ^ 
[c,d] C (0,2) and /i : R — (0, 1) be continuously differentiable. Let {Ti)i>i be a sequence of 
arrival times of a Poisson process with unit arrival time, (Vi)j>i be a sequence of i.i.d. random 
variables with distribution rh^dx) = ^j"^ j~^l[_j^_j4.i[u[j_i,j[(x)(ix on R, and (7t)i>i be a 
sequence of i.i.d. random variables with distribution P{'~fi = 1) = -P(7j = —1) = 1/2. Assume 
finally that the three sequences (r.j)j>i, (Vi)j>i, and {■ji)i>i are independent and define for 
t e R 

(5.20) 

The process Y is h{u)-localisable at all m G R, with = b{u)La(u),h{u) (the well balanced 
linear fractional stable motion). 

Proof 

We apply Theorem 14 . 5 1 wit h m{dx) = dx, r(x) = ^ J2'jLij^'^l-j-j+i[u[j-i,j[{x), f{t, u, x) = 
\t - _ \^\h{u)-i/a{u) g^^^ ^Yie random field 

2-2 °° 

x{t,u) = b{u)cl{:l^\^y/'^^-^ Yl 7.r-^/"(")(|t-F,r(«)-V"W-|T/,l'^w-V"(«))ij_^._^.^^j^j^._^^^. 

X{.,u) is the {a{u) , h{u))-we\\ balanced linear fractional stable motion and it is 
localisable with local form X'^{.,u) = b{u)La(u),hC>j-) El- 



(Csl) The family of functions v —>■ f{t,v,x) is differentiable for all {v,t) in a neigh- 
bourhood of u and almost all a; in ii^. The derivatives of / with respect to u read: 

/:(t, w, x) = {h\w) + ^T^) [(log \t-x\)\t- - (log |x|)|x|'^('")-^/"('")] . 

(Cs2) In [9j, it is shown that, given -u G R, one may choose e > small enough and 
numbers a, 6, with < a < a{w) < b < 2, < h_ < h{w) < < 1 and 
i — I < /i_ < /i+ < 1 — (i — i) such that, for all t and w ixiU := {u — e^u + e) and all 
real x: 

\f{t,w,x)l\fi(t,w,x)\<h{t,x) (5.21) 

where 



Cimax{l,|t-x|''--i/" + \x\^--^'°'] {\x\ < 1 + 2 maxjec; |t|) 

(|x| > 1 + 2 maxtec; \t\) 



for appropriately chosen constants Ci and C2. The conditions on a, 6, /i+ entail that 
supfgf/ \\kiit, ■)\\a,b < oo and (Cs2) hold. 

(Cs3) is obtained with (15.211) for the same reason as in (Cs2). 
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• (Cs4) 



\f{t,w,x)\og{\f{t,w,x)\)r^'"^ < |/(t,^,x)r('") + |/(t,«;,x)log(|/(t,^,a;)|)r("')l{|/(i,.,.)l>e} 

+ \f{t, w, x) log(|/(t, x) I) r(-)l^|^(,,,,)|<i}. 

Since 

\f{t,w,x)\<h{t,x) (5.23) 

one gets 

W, X) l0g(|/(t, W, a;)|)|l{|/(t,^,a;)|>e} < ki{t, X) log(A;i(t, x))l{\fi^t,w,x)\>e} 

< \ki{t,x)\og(k^(t,x))\ 
|/(t,^,a;)log(|/(t,^,x)|)r('")l{i;(,,^,,)l>e} < \h{t,x)\og{k,{t,xW^'"\ 

Fix 7? > such that l<7y<a + |- a/i+ and A > such that ^ < A < 1. 

i/(i,^,x)iog(i/(t,^,x)i)r(-)i{i^(,,^,,)i<i} < 

and thus (Cs4) holds. 

(Cs5) For j large enough (j > j*), 

|/(t,w,a;)log(r(x))|"(^) < i^i|A;i(t, 

+00 

+X2 5]|/(i,«;,x)r(-)(log(i))<^l[ 



Thus 



^, x) log(r(x)) r^-) < K,Mt, x) r(-) + ^4 E T^0^)M-3,-3^-ny^\i-U^) 



j=j* 



To conclude, note that 



j^\^0^)M-3,-3^n-\J-^d^)d^ - 2(log(i))^£^ 



X 

2(a(l + l/6-/i+) -1)- 



(1+1/6-/^+) 

(log(j))' 



a(l+l/6-/i+) 
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Theorem 5.10 (Multistable reverse Ornstein-Uhlenbeck process). Let A > 0, a : R ^ 
[c,d] C (1,2) and 6 : R ^ R"*" be continuously differentiable. Let (rj).j>i be a sequence of 
arrival times of a Poisson process with unit arrival time, (Vi)j>i be a sequence of i.i.d. random 
variables with distribution m(dx) = Y2^=i'^~'^~^'^[-j,-j+iMj-i,j[i^)'^^ ^'^^ (7i)i>i be a 

sequence of i.i.d. random variables with distribution P{'~fi = 1) = -P(7i = —1) = 1/2. Assume 
finally that the three sequences (rj)j>i, (Vi)j>i, and (7i)i>i are independent and define 

+ 00 +00 

yit) = mcl{;''' E E 7.rr^/"(*)2(^-+i)/"We-^(^-*)l[,+..)n([-,,-,+i)u[,-i.))(V^.) (t G R). 
1=1 j=l 

(5.24) 

Then Y is l/a{u)-localisable at any m G R, with local form = b{u)La(u)- 
Proof 

We apply Theorem l4.5l with midx) = dx, r{x) = J2'jLi'^''~^^'^l-j-j+i)uij~i,j){x), f{tjU,x) = 
e~'^*^^'"*^l[i^+oo)(a^) and the random field 

oo 

X{t,u) = biu)Cl{:l-^ E 7.r-^/"(")20-+^)/"(")e-^(^-*)l[,+^)n([-.-.+i)u[,-i.))(^.)- 

X{.,u) is the symmetrical a('u)-reverse Ornstein-Uhlenbeck process and is ^^^-localisable 
with local form X'^{.,u) = b{u)La{u) ESJ. 

• (Csl) The family of functions v f{t,v,x) is differentiable for all {v,t) in a neigh- 
bourhood of u and almost all x in E. The derivatives of / with respect to u vanish. 



(Cs2) 
and 



|/(t,^,x)r("') = e-^"('")(^-*)l[i,+oo[(x) 



sup (|/(t,ti;,x)|"(^Va; < / e-^^^^-'Ux 

1 

< — 
- Ac 

thus (Cs2) holds. 

(Cs3) = so (Cs3) holds. 

(Cs4) 

\f{t,w,x)log{\f{t,w,x)\r^'"'> = X''^'"\x-tr^^^e-''''~-^^^^^ 
as a consequence 

sup |/(t, w, x) log(|/(t, w, x)!)!"^""^ 

R«)e-B(«,e) '- 



+ 00 



dx < / AVe"^""rfM 
Jo 

< +00. 
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(Cs5) 

+00 
i=i 

Fix j* large enough such that for all j > f, l[t,+oo[n{[~j~j+i[u[j-i,j[){x) = 
Then 



/ sup x) log(r(x) dx < 

Jr. weB(u,e) '- ^ 



(j + l)'^log(2)' 



Ac 



+CO „j 

+ (j + l)''log(2)M e-^^("-*) da; 



< i:^^^^^P^ + log(2)V-(e^^-l) £(j + l)V 



-Xcj 

Ac ' "'^'^ 



6 Finite dimensional distributions 

In this section, we compute the finite dimensional distributions of the family of processes 
defined in theorem I4.5[ and compare the results with the ones in [9]. 

Proposition 6.11 With notations as above, let {X{t,u),t,u G R} be as in I 14-1S^ and 
Y{t) = X{t,t). The finite dimensional distributions of the process Y are equal to 

E e = exp -2 / / sm\Y.^,h{t,)-^^^f{t,,t,,x)) dy m{dx) 



E Jo 



/ormGN,0 = GR'",t= 

Proof. Let m G N and write (f)t{0) = E e •'-^ j . We proceed as in [19], proposition 

1.4.2. Let {t/jjjgN be an i.i.d sequence of uniform random variables on (0, 1), independent of 
the sequences {72} and {Vi}, and g{t,u,x) = 6(n)C^^^|"V(a;)"^/"^"^/(t, u, x). For all n G N, 



j=l k=l j=l ^ ^ k=l 

(6.25) 
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m 

The right-hand side of (16.251) converges almost surely to Yl when n tends to infinity 

i=i 

and thus 

(f)t{e) = lim E e 

Set (j/l{9) = E \ e j . This function may be written as: 

\k=l 

All the sequences {7fc}, {Uk}, {Vk} are i.i.d. As a consequence, 

We compute now the expectation using conditioning and independence of the sequences 
{7fc}, {f/fc} and {Vu}. 

E e -1 = E E e |f/i,Vi 



/ m 

E I cos(5^ e,n-V"fe)f/-V"fe)^(^^.^ ^^.^ ^^)) 

- / cos(^^,i/-^/"(*^)^(t„t„V^i))rfi/ 
1- - / E 



The function sin^ is positive and thus, when n tends to +cx3. 



E (^sin^(x: t,, j - E (^sin2(x: I 



To conclude, note that 

\ r Q. 



Comparing with proposition 8.2, Theorems 9.3, 9.4, 9.5 and 9.6 in [9], it is easy to prove 
the following corollary, which shows that the approach based on the series representation and 
the one based on sums over Poisson processes yield essentially the same processes: 
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Corollary 6.12 The linear multistahle multifractional motion, multistahle Levy motion, log- 
fractional multistahle motion and multistahle reverse Ornstein- Uhlenbeck process defined in 
section have the same finite dimensional distributions as the corresponding processes con- 
sidered in 19]. 

7 Numerical experiments 

We display in this section graphs of synthesized paths of some of the processes defined 
above. The idea is just to picture how multistability translates on the behaviour of random 
trajectories, and, in the case of linear multistable multifractional motion, to visualize the 
effect of both a varying H and a varying a, these two parameters corresponding to two 
different notions of irregularity. The synthesis method is described in [8j . Theoretical results 
concerning the convergence of this method will be presented elsewhere. 

The two graphs on the first line of Figure[T]((a) and (b)) display multistable Levy motions, 
with respectively a increasing linearly from 1.02 to 1.98 (shown in (c)) and a a sine function 
ranging in the same interval (shown in (d)). The graph (e) displays an Ornstein- Uhlenbeck 
multistable process with same sine a function. A linear multistable multifractional motion 
with linearly increasing a and H functions is shown in (f). H increases from 0.2 to 0.8 and a 
from 1.41 to 1.98 (these two functions are displayed on the right part of the bottom line). The 
graph in (g) is again a linear multistable multifractional motion, but with linearly increasing 
a and linearly decreasing H. H decreases from 0.8 to 0.2 and a increases from 1.41 to 1.98 
(these two functions are displayed on the left part of the bottom line). Finally, a zoom on the 
second half of the process in (f) is shown, that allows to appreciate how the graph becomes 
smoother as H increases. 

In all these graphs, one clearly sees how the variations of a translates in terms of the 
"intensity" of jumps. Additionally, in the case of linear multistable multifractional motions, 
the interplay between the smoothness governed by H and the jumps tuned by a indicate that 
such processes may prove useful in various applications such as finance or biomedical signal 
modeling. 
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a: Levy multistable motion b: Levy multistable motion 




0.5 1 1.5 2 0.5 1 1.5 2 

xio'* xio"^ 
g; linear multistable multifractional motion zoom on the second half of the process in (f) 




0.5 1 1.5 2 0.5 1 1.5 2 

xlO'* xlO'* 



Figure 1: Paths of multistable processes. First line: Levy multistable motions with sine (a) 
and linear (b) a function. Second line: (c) a function for the process in (a), (d) a function for 
the process in (b). Third line: (e) multistable Ornstein-Uhlenbeck process with a function 
displayed in (c), and (f) linear multistable multifractional motion with linear increasing a and 
H functions. Fourth line: (g) linear multistable multifractional motion with linear increasing 
a function and linear decreasing H function, and zoom on the second part of the process in 
(f). Last line: a and H functions for the process in (g) (left), a and H functions for the 
process in (f) (right). 



36 



